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We discuss the effect of superimposing multiple sources of resonant non-Gaussianity, which 
arise for instance in models of axion inflation. The resulting sum of oscillating shape contribu- 
tions can be used to "Fourier synthesize" different non-oscillating shapes in the bispectrum. 
As an example we reproduce an approximately equilateral shape from the superposition of 
(9(10) oscillatory contributions with resonant shape. This implies a possible degeneracy 
between the equilateral-type non-Gaussianity typical of models with non-canonical kinetic 
terms, such as DBI inflation, and an equilateral-type shape arising from a superposition of 
resonant-type contributions in theories with canonical kinetic terms. The absence of oscilla- 
tions in the 2-point function together with the structure of the resonant Appoint functions, 
imply that detection of equilateral non-Gaussianity at a level greater than the PLANCK 
sensitivity of J'nl ~ 0(5) will rule out a resonant origin. We comment on the questions 
arising from possible embeddings of this idea in a string theory setting. 
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1 Introduction 



Recent years have seen the advent of precise observations of the cosmic microwave back- 
ground radiation (CMB) [Tj, [2J [3] , and the redshift-distance relation for large samples of 
distant type IA supernovae [U [5] , as well as a host of additional increasingly precise measure- 
ments such as baryon acoustic oscillations [6j or the determination of the Hubble parameter 
H by the Hubble Space Telescope key project [Tj. These results are so far in concordance 
with a Universe very close to being spatially flat, as well as with a pattern of coherent 
acoustic oscillations in the early dense plasma which was seeded by an almost scale-invariant 
spectrum of super-horizon size curvature perturbations with Gaussian distribution. This 
structure is a direct consequence of a wide class of models of cosmological inflation driven 
by the potential energy of a single scalar field with canonically normalized kinetic energy in 
(strict) slow-roll. In these models of inflation the spectrum of scalar curvature perturbations 
generated during the quasi-de-Sitter phase is necessarily Gaussian, with 3-point interactions 
generating very small non-Gaussian deviations Jnl = 0(e, if) of order of the slow-roll param- 
eters during inflation [5]. Consequently there are large classes of inflationary models which 
generate sizable levels of non-Gaussianity and are characterized by departures from one or 
more of the defining properties of canonically normalized single-field slow-roll inflation. 

Non-Gaussianity can arise from a non-canonical kinetic term for a single scalar field 
characterized by the presence of higher-derivative terms, from controlled non-permanent 
violations of slow-roll such that the overall inflationary behaviour is preserved, or from the 
presence of several light scalar fields during inflation. A series of results have established 
that these different mechanisms of generating a large 3-point function and its associated non- 
Gaussianity lead to different 'shapes', the distributions of the amount of non-Gaussianity over 
the momentum-conservation dictated triangular domain of the three fluctuation momenta in 
the 3-point function. These shape functions allow for discrimination among many different 
classes of inflationary models, provided one has sufficient measurement accuracy for the 
shape, and one has shown the existence of a complete orthonormal system of shapes of 
non-Gaussianity in the space of inflation models. For two recent reviews see e.g. [U [TU] . 

In this article, we explore the degeneracy, at the level of non-Gaussianity, of different 
classes of inflationary models. Specifically we study degeneracy between models of inflation 
characterized by so-called non-canonical kinetic terms such as DBI inflation [Til H2] and 
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theories which are purely canonical. By non-canonical kinetic terms we mean terms of the 
form X n , where X = |<9 M 0<9 M 0, so that a non-canonical Lagrangian is given by some general 
P(X,(p) [13]. These are thus a subset of the theories described by the Horndeski action, 
which is the most general action for a single field with at most two time derivatives acting 
on it at the level of the equation of motion, i.e. the most general action for a single scalar 
field without ghosts or Ostrogradski instabilities (See e.g. [H] and [15]). 

Theories in this class give rise to an equilateral-type non-Gaussianity inversely propor- 
tional to the reduced speed of sound, c s (see e.g. [16]): 

fz i1 ~ 4 (1-1) 

c s 

where 

c] = (l + 2X^j \ (1.2) 

By contrast, canonical slow- roll single scalar field models give rise to Gaussian spectra 
(with corrections of order the slow-roll parameters) for which all odd-point correlation func- 
tions vanish and all even-point correlations are given in terms of the two-point functions. 
How then could there be any observational degeneracy between a model of non-canonical 
inflation and a canonical single scalar field model of inflation? In fact, single scalar field 
models with canonical kinetic term give rise to (approximately) Gaussian spectra only as 
long as the slow-roll conditions are always satisfied. As argued in [17J, these theories can 
generate observably large non-Gaussianities when there are features in the potential. The 
two classes of such models considered in [T7] are: 1) potentials with steps resulting in a 
localized violation of slow-roll, and 2) potentials with a small modulation which imply a 
large variation rj for small e and rj. Because one of the terms in the three-point function 
is proportional to erj, this can give rise to an appreciable non-Gaussianity with a specific 
shape, termed resonant non-Gaussianity [IS], which is orthogonal to the local and equilat- 
eral shapes. Resonant non- Gaussianity comes with the additional property of allowing for 
very efficient calculability of the iV-point functions up to iV = 10 ... 20 |19j . 

However, as we shall show, a closely related model, in which a sum of such modulations 
in the potential is present, leads to a very different non-Gaussianity. Its shape is of the 
equilateral type in the sense that its momentum average has a large overlap with the standard 
equilateral one, despite the lack of non-canonical kinetic terms. Yet, periodic scale dependent 
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features remain. Thus observation of equilateral non-Gaussianity might not necessarily imply 
that non-canonical kinetic terms were present or relevant in the underlying model of inflation. 
This degeneracy at the level of the three-point function leads us to conclude that it may 
not be as easy as previously hoped to distinguish amongst single-field models of inflation. 
However, we then clarify that the non-observation of oscillating contributions to the 2-point 
function, and the structure of the resonant iV-point functions place tight constraints on the 
maximum amount of resonantly generated non-Gaussianity with equilateral characteristics. 
In particular, this implies that detection of equilateral non-Gaussianity at a level greater 
than the PLANCK sensitivity of Jnl ~ 0(5) will rule out a resonant origin [20]. In case 
a detection of non-Gaussianity eludes PLANCK, future 21cm observations may provide a 
detection capability of fai > 0(0.01) down to the slow-roll level itself due to the extremely 
large number of modes (0(1O 16 )) available [21, 22J. For equilateral-type non-Gaussianities of 
this magnitude there is then a degeneracy between a possible non-canonical and a resonant 
origin. 



2 A sum of resonant bispectra 

2.1 Potential, solution and slow-roll parameters 

From an effective field theory point of view, we can begin with a model of a single scalar 
field with canonical kinetic term and modulated potential 

V{cj>) = Vo{4>) + cos (^jf) , (2-3) 

where we have generalized the modulated potential in [TH] to the case where the modulation 
is a series of terms with varying phases. For suitable coefficients and values of /j, the sum 
remains a small perturbation on the potential Vq(4>) (see Figure [TJ. We will show that slow- 
roll inflation for this theory is supported for sufficiently large /j, while a large variation in 
the slow-roll parameters is possible. In Sections 2.2| and 2.3 we give the power spectrum 



and bispectrum respectively, finding that the bispectrum is given by a series of the resonant 
bispectra found for a singly modulated potential in [18] . We shall follow closely the analysis 
in [18]. The details of our calculations are given in Appendix [Aj 



The equation of motion for the inflaton with potential eq. (2.3) is 
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Ji V Ji 



(2.4) 



As in [TS] we treat the oscillations as a perturbation and write <fi = 0o + <\>\ + ••• where 
is the solution in the absence of modulations, given by 

A Wo) 



(2.5) 



The equation of motion for 1; written in terms of derivatives with respect to O (denoted 
by a prime) is 

3 3?7y \ 3 ^ M ■ ( 0o + c 



where Vo = Vo(0o)> an d ey and 77^0 are the canonical slow-roll parameters evaluated for the 
potential V . Note that e' Vo = y/2e Vo (r] Vo — 2e Vo ). The solution is a generalization of that 
found in [18]; to first order in slow-roll parameters and assuming /j <C y/2e±\/i, it is given 
by 



w<) =-dw?^" h (^i' (2 - 7) 



where the subscript * indicates that e = ey o (0*) is evaluated at horizon exit. This solution 
works because 0i = J^0i, where each 0j satisfies the EOM for a specific /j. 0o(^) is time 
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dependent but, for the period of interest, in which modes observed today exit the horizon, 
we can approximate it as being given by the value at horizon crossing, i.e. 0*, except in the 
oscillating term. In other words we have assumed that the dominant variation in <f>\ is due 
to the modulating terms. Defining 

we have that the full solution is well approximated by 

*(f) = 0o(t)-£^sin(^). ( 2 - 9 ) 

The slow-roll parameters e = — and 5 = can ^ e calculated for the potential 



eq. (2.3) and are found to be: 

2 



Vn 



1) 



e ,-3v^E^ cos (^) 



e = e + ei + 0((b* 
8 = e* - 77* - v / 2e^0i + ~ £*)4>u 



(2.10) 



5 = 5 o + 5i + 0((&-) 2 - 

The SR parameter of interest for observation is ^ where 5q = — 77* is the value of 
5 coming from the unmodulated potential Vq. A large jj is what is responsible for a large 
resonant non-Gaussianity in [18]. Here, we find 



Note that : ^¥ ± is large for all i for the solution eq. (|2.9|). 



2.2 The power spectrum 



As in [T5], the Mukhanov-Sasaki equation for the mode function TZ k of the curvature per- 
turbation in our case is given by (see the Appendix for details) 



d 2 TZ k 2(1 + Si(x))dK 



k 



dx 2 x dx 
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+ n k = o, (2.12) 



where x = —kr for r the conformal time. We have neglected eo and Sq and are working only 
to leading order in e*. The form of this equation is the same as in 1181 . but Si is now defined 



in eq. (2.10) and given by 



, 1 = -3£«f.in(*^) 



(2.13) 



We will find a sum of integrals each encountering resonance at a value £ res ,i- As in the 
case of a single modulation, the effect of 8% is negligible both at early times x 3> x reS)i and 
at late times x <C x res i. We therefore look for a solution of the form 



TZk(x) 



K 



(o) 
fc,0 



/ , / -x 3/2 H ( M(x) - cl ] (x)i *lZx 



^Hi%(x) 



(2.14) 



where TZ^ is the value of TZk{x) outside the horizon in the absence of modulations, H^\x) 
are Hankel functions and c^~\x) vanishes at early times. Then the Mukhanov-Sasaki equa- 
tion gives an equation governing the time evolution of cjj. (x) which is given by (for fa <C y/2e^ 
for all i) 



d 

dx 



-2ix_^_ (-)/ \ 

6 dx Ck [X) 



-2i 



Si(x) 



x 



(2.15) 



To solve this we need Si as a function of the conformal time r. In terms of r the 
background equation of motion is given by 



(2.16) 



dln(-r) 

with solution 0o( r ) = 0* + \/2e*ln ^. This can be written as 0o( a; ) — 0fe + \/2e*lnx, where 
0fc = 0* — \/2e*ln jr- is the value of the field when the mode with comoving momentum k 
exits the horizon. Then, writing c 



(-)/ 



c i A 3 ')) we nee d to solve 
2^7 



72T cos( 7; + 7; + ^ 



ma;). 



(2.17) 



Thus each c^^x) behaves like c ( k '(x) in [18] but with b* = b\, f = /, and resonance at 



,(-) 



different x values x 



res, « 



For each term we can use the stationary phase approximation to find 

fi ^ 



,(-) 



U k (x) 



3 i(^) 
e /i , 



(2.18) 



! T 3/2fr(l), 



)-E 



3i6*7r 



i — i( 

e 



x^Hf/ 2 {x)].{2.lV) 



Outside the horizon, k/aH <C 1 i.e. x = —kr = k/aH <C 1 (using the definition of 
conformal time and the fact that during inflation H is constant). Then we find 



and 



with 



n k (x) = 



i 



fan ( 0fc + Ci 
cos 



2y'2r 1 V ft 



: sin 



n 



+ 0{x i 



1/2 



COS 



*>k + c, 

/i 



which is just the generalization of (2.30) in [TS] one might have expected. 



(2.20) 



(2.21) 
(2.22) 



2.3 The bispectrum 

Next we ask if this generalization carries through to the bispectrum as calculated in [18J. As 
in the case of a single modulation, the leading contribution comes from 



H T (t) D - / d 6 xa 6 {t)e{t)5{t)K l {x,t)K{K,t)i 



(2.23) 



where to linear order in b* we can use the approximations e ~ e*, 5 ~ 5\ and use the 



unperturbed mode function^] 



n k {x) = n^-x^H^ 



X). 



(2.24) 



The analysis of [18] therefore carries over to our case, with the three-point function given 



by 



< ^(k!, t)ft(k a , t)K(k 3 ,t) >= (2vr) 7 A^-^r2^( k i + k 2 + ka! (J ^ h) 



k 2 k 2 k 2 
^1 ^2 3 



(2.25) 



for 



hk 2 k 3 8 



dX-^-e 
H 



1 ^ ki i K{k\ + kl + k\) 
X ^~k~ + X 2 hk 2 h 



+ c.c, (2.26) 



1 To see why we can use the unperturbed mode functions, see Appendix A. 2 
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where K = k\ + k 2 + k 3 , X = —Kr, and <5i is given by eq. (2.11). Equivalently this can be 
written 



g(fci,fc 2 ,fc 3 ) _ 1 
where 



E / dx 



St 
H 



i -iX 



~* ~~ X ^ h; + X* k x k 2 k 3 



+ c.c. (2.27) 



IT 



* 3b* cos 



^>0 + Cj 

fi 



—-^35* cos 

Ji 



> K + y/2e7lnX + Q 



(2.28) 



2r, 



Now, for each integral in the sum the biggest contribution comes from X resi — 

' Ji 

can be seen by defining 



This 



1-KA 



h 

36*x/2^ 
2/i 



dXe" iX cos 



'A" 



+ v/2eJnX + c s 



(2.29) 



The phase of the first term can be zero while the phase of the second cannot, which 



means that the first term will dominate, with its dominant contribution at X ms i — — . 

' Ji 

This confirms that the main contribution to the bispectrum arises when the modes are still 
deep in the horizon. In this regime, we can approximate 1/X and 1/X 2 by 1/X res ^ and 
1/X r 2 esi in eq. (2.27) and we have 



Q(h, k 2 , k 3 
hk 2 k 3 



T T 1 V- h 1 K{kl + kl + kl) ^ 

^res,i „ , Km ^res i K\K 2 K 3 



Further we can evaluate Xka in the stationary phase approximation to get 



36*V2vr ( JTe 



-KA 



fi 



3/2 



;! (fa+e,)//, e -f e -(f-f l"f )^ 



(2.30) 



(2.31) 



2e*lnf , this gives us the shape of the non-Gaussianity as a sum of 



Using 4> K = 0* - 

the resonant non-Gaussianity shapes in [18] (up to an overall phase) j^] 



2 Note that each term in the sum over i in eq. ( |2.32| actually has a different phase contribution which 
depends on the fi, see eq. (2.31 ). However, we can absorb these terms into the which we are free to choose, 



such that eq. (2.32) is correct. 
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Gres(ki, k 2 , k 3 



E 



3V27T&* ( JTe 



fi 



3/2 



sin 



')_K+_Ci 
fi 



fi ke 
^ — cos 



2c* 

£=f=m 



fi y K(k 2 l + k 2 2 + k 2 3 ) . + q 

1 " — sin 1 



E 



'2e 



3/2 



sm 



fi 



^ln- + ^ | 

/i ^* fi 



fi sr^ k, 



E 



^2e* jr^ "V. 



COS 



In — + -i I -f ... 



/i ^* fi 



/i 



(2.32) 



As a consistency check of our calculation, we have checked that this sum of resonant 
bispectra satisfies the squeezed limit consistency relation [SJ [23] • The details are presented 



in Appendix |A.3 



3 Equilateral features from summed resonant non-Gaussianity 

We now want to discuss what kind of non-Gaussianities could have a sizable overlap with 
summed resonant non-Gaussianity. One could also turn this question around and ask what 
values the parameters b*, fi and q have to take in order to generate a degeneracy. From 
an effective field theory point of view it is perfectly fine to choose ad hoc values for these 
parameters as long as they are not in conflict with observations such as the power spectrum 
and fulfill certain consistency conditions. For instance, the frequencies should certainly not 
be super Planckian, i.e. fi < 1, and monotonicity of the inflaton potential requires b* < 1. 
A possible stringy origin of these parameters is discussed below in Section |4j 

In bounds on oscillating features in the power spectrum were given for a quadratic 

and a linear inflaton potential respectively. The bound that was found in both works is 

10" 5 

b*fi < -= , (3.33) 



for a single modulation in the potential. We assume that this bound is valid for multiple 



modulating terms in the potential. Regarding the applicability of the bound eq. (3.33) to 



general potentials Vq (</>), see the discussion in Section 3.4 
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Let us briefly discuss the general form of the bispectrum. Due to the appearance of 



the delta function in eq. (2.25), a momentum configuration is completely characterized by 
the absolute values of the three momenta ki, k 2 and k%. Furthermore, for a scale-invariant 
spectrum this reduces to two variables. Then, one usually considers the bispectrum as a 
function of the two rescaled momenta X2 = k 2 /ki and x$ — k^/k\. A region that includes 
only inequivalent momentum configurations is given by 1 — X2 < £3 < Xi- 
Let us switch to the variables x± = X2 ± #3, with 1 < x + < 2 and < x_ < 1. Note that the 



resonant non-Gaussianity eq. (2.32) is to first order in f i /y/2e+ only a function of x + and k\ 
but not of X- since 



sin 



— In- 

fi k* 



sm 



fi 



(y + In fa/k* 



(3.34) 



having defined y = ln(l + x + ). Note that summed resonant non-Gaussianities are therefore 



not scale invariant, as is clear from in the explicit dependence on k\ in eq. (3.34). We will 



discuss the issue of scale dependence in Section 3.2 



Furthermore, eq. (3.34) implies that as far as other types of non-Gaussianities are con- 
cerned, we can only expect degeneracies with the summed resonant type if they are pre- 
dominantly a function of x + . We will show in the next section, Section 3.1| that this is 
primarily the case for equilateral non-Gaussianity, typically arising in non-canonical models 
of inflation. 

To measure the degree of degeneracy between different kinds of non-Gaussianities we fol- 
low [2S]- The bispectrum can always be parametrized in the form g ^^ 2 fc '^ 3 - > ~ fNzS(ki, k 2 , k 3 ), 
with 'amplitude' f^L and shape function S(ki,k 2 ,k3). In general, the shape refers to the 
dependence of S(k\, k 2 , k%) on the momentum ratios k 2 jk\ and ks/ki when the overall mo- 
mentum scale K is fixed, while the dependence of S on K when the momenta ki are fixed 
gives the running of the bispectrum. 

Now the cosine of two shapes is defined via the normalized scalar product 

F{3, S') 



C(S, S') 



^/F(S,S)F(S',S'y 



where 



with 



F{S, S') 



dV 



dV 

~K 



K 

dki dk 2 dk s 
h + k 2 + k 3 
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S(h, k 2 , k 3 )S'(ki, k 2 , h) , 



1 



-kdk da df3 , 



(3.35) 

(3.36) 
(3.37) 



where we have switched to the variables k, a and (3 defined according to 

k = -(k 1 + k 2 + k 3 ), fci = fc(l-0), 

k k 
k 2 = -(1 + a + p), h = -(1-a + p), 



(3.38) 



with the integration boundaries 

a e [-1 + 0, 1 - P] , Pe[0, 1] and k e [k 

mini ">max 

] . (3.39) 

3.1 Equilateral and local shapes for x _ — > 

Out of the many known types of non-GaussianitiesJ^] we discuss the following two represen- 
tative types: The equilateral type 

Gequil(ki, k 2 ,k 3 ) _ „equil 



hk 2 k 3 

with 



fNL S equil(h, k 2 , h) , (3.40) 



Q n u i\ {k 1 + k 2 -k 3 )(k 1 + k 3 -k 2 )(k 3 + k 2 -k l ) 

Sequilih, k 2 , k 3 ) = — — , (3.41) 

K\K 2 K 3 

is characteristic of non-canonical inflation. 
The local type is given by 

Giocai(ki,k 2 ,k 3 ) _ local 



k\k 2 k 3 

with 



f%fS local (h,k 2} k 3 ), (3.42) 



k x k 2 k 3 

and is dominant in multi-field inflation for instance. 

In the limit x_ — > 0, the shape functions are given by 

-ix--+0/__ \ 4(x + — 1) 



1*3 I 1.3 I h.3 

S local (h,k 2 ,k 3 ) = kl Tf 2 , 3 , (3.43) 



S X eq - u 7M = V \ \ (3-44) 
Sf -J°(x + ) = ^. (3.45) 



Now, one can evaluate the cosine eq. (3.35) to find the overlap of S equ u and S* ^{x^ 



C(S eouih S x e -J°)= 0.93. (3.46) 



5 For an overview see e.g. [9]. 
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Hence, the equilateral shape is well approximated by its x_ — > limit. For the local shape the 
overlap is much smaller. Si oca i diverges for squeezed momentum configurations which makes 



it necessary to regulate the integrals that enter the cosine eq. (3.35). We only give an upper 



bound C (S i oca i, Sf oc ^°) < 0.7, which was obtained by cutting the integration boundaries 



eq. (3.39) as follows: 



p e [A, 1 - A] with A = 5 • 1(T 5 . (3.47) 

The different qualitative levels of agreement of the equilateral and local shape are also 
visualized in Figure [2j Notice that the shape function S non - can from non-canonical inflation 
shows a slightly more complicated momentum dependence than the equilateral shape func- 



tion eq. (3.41), see [16J. We find C(S non _ can , S^ /0n _ can ) > 0.93, so the approximation by the 



,x_ — > limit is even better in the case of S. 



non-can- 



3.2 Scale invariance vs. scale dependence 

Let us now discuss how well the scale-invariant equilateral shape can be approximated by 
a scale-dependent shape, such as summed resonant non-Gaussianity. It is obvious that the 
overlap cannot be made arbitrarily large; however, we will show in the following that the 
overlap can still be considerable. 

Let S v ^ uil (y) be the periodic generalization of 

S: q - U u°(y) = 4 (|^ , V G [In 2, In 3] , (3.48) 

to y e M, i.e. S^Z u {y + Ay) = S^ uil {y) with Ay = In 3 — In 2. This definition is motivated by 
the fact that we are going to Fourier expand on the interval y e [In 2, In 3] and this expansion 
will itself be periodic with period Ay. 

Furthermore, let us consider the scale-dependent shape S^Z u (y+hi ki/k±) = S^ u (\n K / k± 
which is a shape that can be approximated by the shape of some combination of resonant 



non-Gaussianities, since they have the same functional dependence, see eq. (3.34). The 
different shapes in (t/,ln£;i//c*) space are shown in Figure |H 



The calculation of the cosine eq. (3.35 ) between the shapes S^ u (\n 2k /k+) and S equ u(a, (3) 
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Figure 2: The equilateral and local shape functions (left) and their approximations x_ — > 
(right). 

simplifies due to the reduced dependencies on the integration variables k, a and (3 to 



r< ( q qP er \ 

^ \Jequih J equil) 



J da d/3 S equil (a,/3) 



2j k k :: dkks Zu^2k/k 



e 1 uil > / \ 1/2 / , „ \ 1/2 

2 / da dp S^ 2 (a, /?)) - **J dk k S^ 2 (ln ^/K 

(3.49) 



We find that eq. (3.49) is to very good approximation independent of the values of k min and 
k maxi if k min < k* and k max > k±. This is due to the periodicity of Sl^ uil {hi2k/K), yielding 



the second fraction on the RHS of eq. (3.49) essentially independent of these values. For 
kmax = 10 4 fc* and k min = 10~ 4 £;*, we find 

C(S equih S%: u )= 0.83. (3.50) 
Let us also mention that the overlap would be greater than 83%, if we had considered 
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Sequil* "*°(y) 



Se q uii per (y+ln(fci/fc,)) 




Figure 3: The scale-invariant shape S x equil (y) (left) and scale-dependent shape S^ u (y + 
lnki/k*) (right) for y G [In2,ln3] and ln&i/A;* G [0, Ay] 



not the equilateral shape function eq. (3.41), but the exact shape function [16] that arises in 
non-canonical single field models of inflation. This is due to the fact that the overlap between 



S x -^° and S is larger than for the equilateral shape function eq. (|3.41|), as discussed at the 



end of Section 13.11 



3.3 Fourier analysis 

Having found that the equilateral shape function has a non-negligible degeneracy with the one 
dimensional function S^Z a (y + hiki/k*), the question of degeneracy with summed resonant 
non-Gaussianities can be expressed in terms of a Fourier analysis. Let us set h\ = k* for 
conciseness in the remainder of this section; however the following analysis is valid for all 
values of k\. 

Let us define 



= -^4— l ( y and F * = • ( 3 - 51 ) 



Then, summed resonant non-Gaussianity, eq. (2.32), can, to first order in fi/ y/2e+, be written 
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in the forrrQ 

Qre 



2N 



N 



hk 2 k 



(3.52) 



i=i 



i=l 



where in the last equality of eq. (3.52), we have chosen the phases Cj such that there appears 
a sine and a cosine for each frequency Fj. 

The Fourier expansion of a generic function f(y) on the interval y G [a, a + b] is given as 

A? 



f(y) 



2 



E 



cos 



2vu 



«Af+i sin 



2iri 



-y 



(3.53) 



with 



UN+i 



a+b 



f(y) cos 



a+6 



f(y) sin 



27ri 
IT 

~b~ 



y J dy for z > , 
y I dy for z > 1 . 



(3.54) 



Note that eq. (3.52) is of the Fourier form, eq. (3.53), if the constant term 



U 



f{y)dy 



(3.55) 



is zero. Since S^ u (y) is monotonically increasing and has a zero at In 2, the numbers a and 
b can always be chosen such that u = 0. Note that we can always redefine S^ a (y) on the 
(non-physical) interval [In 2 — £,ln2] for £ 1 such that the integral eq. (3.55) has a large 
negative contribution on this interval, i.e. 



u 



a+b 



f{y)dy 



In 2 



In 2-5 



In 3 



f(y)dy+ / f(y)dy 



(3.56) 



In 2 



where the first term on the RHS of eq. (3.56) is negative and has the same absolute value 
as the positive second term such that Uq = 0. For instance one could use f(y) = —\fo\ with 
|/o] >1 fory E [In2-£,hi2]. Hence, b = In 3 -In 2. 

In Figures |1] and |5j we show the Fourier expansion of S^ u (y) for iV = 5 and N — 10. 
The coefficients Ui are given in table [T] for N = 5. 

Due to the suppression by fi/y /, 2e ir <C 1, the second order contribution to Qres/(k\k2kj) induces a 
non-negligible correction only in the squeezed momentum configurations which we ignore in the following. 
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0.5 




Figure 4: Fourier expansion of S^uty) for iV = 5 (left) and N = 10 (right). 



Ml 


u 2 


u 3 


U4 


u 5 


Uq 


M 7 


u 8 


M 9 


Uio 


0.41 


0.01 


-0.17 


0.04 


0.09 


-0.20 


0.24 


-0.03 


-0.12 


0.05 



Table 1: Fourier expansion coefficients Ui, eq. (3.54), of S^ uil {y) for iV = 5. 



Let us discuss how the Fourier expansion parameters a, b and Ui of eq. (3.54) relate to 



the parameters Bi and Fi of eq.'s (3.51) and (3.52). Clearly, if we want to impose 



Geq r uil(y) ~ Qres{y) 



we have to set 



B^/X'«,<0.41f 



NL 



and 



2tu 
~b~ 



15.5z 



(3.57) 



(3.58) 



where we have used that the maximal Ui is Mi = 0.41 for N = 5. In the derivation of the 



bispectrum in Section |2.3| we assumed /, <C v^^*; i- e - ^ 1- We see from eq. ( |3.58[ ) that 

J NL 



this constraint is satisfied with growing i. Now, solving for b* and fi in eq. (3.51) yields 

pequil 



8 £ 



equil 
NL 



11; 



< 



(max 
1 NL 



and 



1 3v/2¥ i^ 2 " f. 
where we have once again used max{wj} = 0.41 and / 



fi 



max 
NL 



/2e* 
Fi 

140. 



0.06- 



x 2e, 



(3.59) 



There are two things we can conclude from eq. (|3.59|). First of all, frequencies are 

d, we see that for f e N q f > f™ x 
140 is the maximal equilateral 



necessarily sub-Planckian as is required for consistency. Second, we see that for f e ^ 1 > f^l 



the monotonicity condition b\ < 1 is violated. Hence / 



max 
NL 
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non-Gaussianity that could be matched by summed resonant non-Gaussianity if there were 
no constraints from the power spectrum. The Jnl that is consistent with the data is much 
smaller than 140, as we will discuss in the following section. 



3.4 Constraints from the power spectrum 



S per (U 2 ^ 3 ) 



1.0 




Figure 5: Fourier expansion of S^ u (y) for iV = 5 (left) and N = 10 (right). We have set 
ki = h. 



Let us discuss the constraint on the product b*fi from the power spectrum, i.e. eq. (3.33 ) 



Using eq. (3.59), we find 

btfi< 



tW f equil V 



NL 
,,5/2 



< 5 ■ 10" 4 /, 



4 fequil /K 
NL V Ze * j 



(3.60) 



where the product on the far right is the maximum value of b*fi. Since the maximum value 



of b*fi has to be smaller than the upper bound given in eq. (3.33) this implies 



e+nf 1 <io- 2 



(3.61) 



Given the values for b* and obtained for the Fourier expansion, we can check the 



behaviour of eq. (2.22). We find that to leading order it behaves like a single oscillation with 



fi given by the lowest frequency fx. This means that the two-point function bounds found 
for the singly modulated potential will still apply in our case. 
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Note, that eq. (3.33) was derived for large-field models (linear axion monodromy infla- 
tion [25J and a quadratic potential [24]), and thus should not have validity for, say, e* < 10~ 3 . 
Therefore, from this result we can argue at most for a resonantly generated f^ U < 0(1). 
More generally, an analysis of inflation in presence of periodically broken shift symmetry 
which generates resonant non-Gaussianity showed [20] that the structure of the resonant 
A-point functions displays a hierarchical suppression with increasing N. Hence, the non- 
observation of oscillating contributions to the 2-point function places tight constraints on the 
maximum amount of resonantly generated non-Gaussianity to no greater than the PLANCK 
sensitivity, that is /jvl ~ C(5). 

4 Possible avenues for embedding: Axion monodromy inflation 

We have seen that a potential modulated by a sum of small oscillating terms as given in 

v(<f>) = VoW + cos i^jf) > ( 4 - 62 ) 



eq. (2.3) 



will lead to a bispectrum given by a sum of resonant bispectrum terms eq. (2.32). This 
follows as long as /j <C i/2e^ for each i. Furthermore, a Fourier series of these resonant 
bispectra leads to a summed resonant bispectrum which closely approximates the equilateral 
bispectrum. Taking this series amounts to a choice of parameters {/i,Cj}, consistent with 
the bounds above, for a suitable number of terms in the bispectrum or equivalently in the 



modulated potential eq. (2.3) 



It remains to ask how such a sum of sinusoidal corrections to the potential could arise. In 
this section we argue for one possible source for these corrections in the context of axion mon- 
odromy inflation [271 ES] , where the resonant non-Gaussianity shape was first observed [18] . 
The oscillatory correction considered in these papers arises from nonperturbative effects such 
as instantons, and is quite general in large field models. We show that these effects can also 



give rise to a sum of periodic potentials such as that in eq. (4.62) and studied here. 



4.1 Axion monodromy inflation 



Up to this point we have taken the potential eq. (4.62) as the starting point, in keeping 



with a low-energy effective field theory point of view of inflation in which we are agnostic as 
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to the UV origins of the theory. One possible avenue for a realization of multiple resonant 
non-Gaussianity is given by axion monodromy inflation [271125] . Axion monodromy inflation 
is a string theoretic realization of large field inflation, where the approximate shift symmetry 
of the axion protects the required flatness of the potential from higher-dimensional operators 
which become relevant when the field range is super-Planckian. Axions in string theory arise 
either as the Hodge duals of 2-form fields in 4 dimensions or as zero modes of p-forms on 
the compact manifold [28], and inherit their periodicity from these fields. Normalization of 
the axion kinetic term results in a coupling to other fields inversely proportional to f a , the 
axion decay constant. In general string theoretic constructions this decay constant is larger 
than the range allowed by cosmological bounds (to avoid excessive production of axionic 
dark matter which could overdose the universe) [28]. However, despite this, the field range 
spanned within a single period is generally sub-Planckian in string theory [29J. 

This limitation (for the purposes of large field inflation) is circumvented in axion mon- 
odromy inflation, where a large field range for the inflaton axion field c = 2ir J 2 Ci is possible 
for instance when the two-cycle £ on which it is supported is wrapped by a space-filling 
5-brane which breaks the shift symmetry c — > c + 2ir in the potential, and gives rise to 
monodromy in the potential energy: the potential is no longer a periodic function of the 
axion, and increases without bound as the axion VEV increases. This follows from the DBI 
action, which gives a potential 

Vb(c) ~ Vi 4 + c 2 , (4.63) 

where l\fa' is the size of the wrapped cycle. Generalizations of this effect arise generically 
in the presence of fluxes. 

The initial axion VEV is large, leading to a linear unmodulated potential Vo(0), which 
means that = the value of the inflaton at horizon exit, and b* = b is independent 
of 0*. Inflation ends when at small values of the axion VEV, Vo(0) is no longer linear, the 
axion oscillates about its minimum, and the inflationary energy is diverted into any string 
modes coupled to the axion. 

Many of the consistency conditions (such as e.g. the absence of large back reaction from 
the axion-induced 3-brane charge building up on the 5-branes) for axion monodromy inflation 
were already discussed in [27J |25]. Some remaining questions may consist of e.g. finding an 
explicit realization on a concrete compact 3-fold in type IIB (despite the generality of the 
mechanism which in the context of fluxes [30] should lead to many explicit avenues for 
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construction), or details of the 3-brane charge back reaction and the mass scale of some of 
the KK modes in particular setups. For a more detailed discussion of some of these and 
related points see e.g. [251 EI]- 

4.2 Periodic corrections to the potential 

In fI7\ |23] , working in an 03/07 orientifold of Type IIB, the axion found to be suitable 
for this model of inflation was that arising from the zero mode of the RR two-form Ci on 
a 2-cycle S in the internal manifold, denoted c = 2n J E Ci. The co homology groups 
split into two eigenstates under the action of the orientifold in a flux compactification [32] 
(see also [33J) 



where the ± subscripts refer to even/odd behaviour under the orientifold action. Because 
C*2 is odd under the orientifold projection, only the mode coming from wrapping a 2-cycle 
which is also odd under the orientifold will survive, i.e. we take S G . Such an odd 
cycle can be written as v~ = v 1 — v 2 , where v 1 and v 2 are two-cycles in the CY which are 
mapped into each other by the orientifold action. The + combination is then obviously an 
even two-cycle. 

Generally, one expects non-perturbative effects to lead to small periodic corrections to 
the potential for the axion, leading to the modulated potential which gives rise to resonant- 
type non-Gaussianity [251 HE]- These corrections arise from instanton corrections to the 
action, specifically worldsheet instanton contributions to the Kahler potential. Instantons 
arise when the (p+l)-dimensional world volume of a Euclidean p-brane is wrapped over 
a (p+1) cycle in the internal manifold. Their corrections are exponentially suppressed by 
the size of the wrapped cycles. Space-time instantons, which correct the superpotential, are 
localized in R A , while worldsheet instantons, which correct the Kahler potential, are localized 
in 2 dimensions. These are due to Euclidean Dl-branes and their SX(2, Z) images (i.e. F- 
and D-strings as well as their bound states, (p, q) strings) wrapping 2-cycles v + threaded by 
the axionic C2 forms. In [25] this correction was taken to be of the form 



for the c axion, where Ve is the Einstein-frame volume of the CY manifold, and v + is the 
even cycle dual to that wrapped by C2 and on which the instanton is supported, r is the 



(4.64) 




(4.65) 



21 



axion-dilaton t = Cq + ie ^ and the G a are complex scalar fields composed of the RR and 
NSNS axionic fields c and b as [M] 

G a = c a_ rb a^ (4 gg) 

where a runs over the number of axions. 



The expression in eq. (4.65) captures the important sinusoidal features of the instanton 
corrections to the Kahler potential. The detailed form of one such series of leading (in the 
sense of exponentially unsuppressed) corrections in the 6-axions was found in [53] . It is given 
by 



g(r, f , G\ G a ) = — V n p V ) cos((n + rm) — ) - mk a G a ), (4.67) 

4*-^ z — ' n + rnr r t — t 

P (m,n) 

where (3 are curves in the negative eigenspace spanned by the basis {w a }, the rip are 
the Gopakumar-Vafa invariants for these curves and K a = jpOJ a - The sum over (m, n) 
corresponds to a summation over SL(2, Z) images of worldsheet instantons, which arise from 
F-strings wrapping 2-cycles. This ensures inclusion of Euclidean Dl-branes and bound states, 
(p, g)-strings, in the instanton corrections. Generically then, at large \bi\ this correction is a 
function of all the frj-axions present. 

On general grounds then (see also the short discussion in [35] on p. 9) we expect the 
presence of a qualitatively analogous series of corrections which depend on the c l -axions at 
large |c*| and which are exponentially suppressed in the volumes of the curves in H^- As these 
terms will again contain a summation over all images under 5X(2,Z), the individual terms 
will not be ~ cos(c*) for a given d but instead should have a dependence ~ cos(a(n, m, r)jC l ). 

We may thus proceed by assuming the corrected Kahler potential to be of the form 

V e - e~^3^ v4 njm (r, f ) cos [a(n, m, r, f)c - ip(m, n, r, f , b)] > , (4.68) 

n,m J 

for the case of one c-axion. The moduli potential often stabilizes b due to its appearance 
in the Kahler potential already at the leading level [3U [27). Note, that this fixes b not 
necessarily at zero which can give rise to a non-zero phase ip(m, n, r, f, b). 

This correction to the Kahler potential results in a proportional correction to the scalar 
potential of the theory, as well as corrections to the kinetic terms of r and G a which are 
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exponentially suppressed in the curve volume v + . Hence, normalization of the axion kinetic 



term, as e.g. in (5.19) of [25], will yield a correction to the potential of the form of eq. (4.62) 

SV ~ e-^+^Y^AiCosictij -ipi), (4.69) 



where we renamed the corresponding coefficients in eq. (4.68) as effective parameters A4, o^. 

For a very crude estimate of the number of terms one may expect to contribute to SV 
we note that already values of n, m smaller than 3 give O(10) possible terms in the SL(2, Z) 
sum, each of which can have contributions from the sum over rational curves in H^- This 
could easily give the O(20) terms required in the sum of resonant bispectra to reproduce an 
equilateral-type non-Gaussian signal. 

Thus we have argued it to be plausible that a series of sinusoidal corrections with suitably 
varying frequencies may be achievable in connection with an axion monodromy generated 
potential, using the instanton corrections already considered at leading order in previous 
work. While the necessary and much deeper analysis of such constructions is left to the fu- 
ture, such a setup could form the starting point for a concrete model realizing the mechanism 
described in the preceding sections. 



5 Conclusions 



We have seen that a potential modulated by a sum of small oscillating terms will lead to a bis- 
pectrum given by a sum of resonant bispectrum terms which can sum to a bispectrum which 
closely approximates the equilateral bispectrum. This has potentially relevant implications 
for the aim of differentiating between inflationary models using data on non-Gaussianities. 
This potential differentiating power rests on (a) non-Gaussianity being one of the few possi- 
ble consequences of inflation that is not generic, (b) different classes of inflationary theories 
giving rise to different dominant bispectrum shapes, and (c) new data allowing us to set 
bounds on the amplitude of different shapes in the bispectrum being expected in the next 
few years [36] . 

In broad terms, single scalar field inflationary models with canonical kinetic terms have 
non-Gaussian signal of order the slow-roll parameters, so that they would be ruled out by 
any measurable detection of non-Gaussianity. P(X, <fi) theories (or more generally, Horn- 
deski theories) with non-canonical kinetic terms of a single inflaton field give rise to pre- 
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dominantly equilateral-type non-Gaussianity [TB], while multi-field inflation gives rise to 
predominantly local-type non-Gaussianity and folded-type non-Gaussianity could indicate a 
non-Bunch Davies vacuum (For reviews see [9], [10] and references therein). Orthogonal-type 
non-Gaussianity, found in [37] , can be the dominant type in multi- field DBI galileon inflation 
with an induced gravity term [38 j [39]. Different contributions to the non-Gaussianity can 
also appear in combination, as in [40J. 

Here we have focussed on the possible degeneracy of equilateral-type non-Gaussianity, 
associated with single field theories with non-canonical kinetic terms, and summed resonant 
non-Gaussianity in a theory with canonical kinetic terms, arising from a potential modulated 
by a series of sinusoidal terms. It has been noted that within the class of non-canonical the- 
ories described by the Horndeski action, it may be difficult to differentiate between different 
models because some linear combination of the same shapes is always present [4"T| l4*2l [14] . 
The best way to differentiate between these models is therefore via conditions on the relative 
amplitudes of these terms which will differ for different models |43j . Despite this, one might 
hope to be able to identify the presence of non-canonical kinetic terms, via the dependence 
of the equilateral signal on the reduced sound speed ft? L ~ \. 

However, there is considerable degeneracy between observables even at this level. A 
dominantly equilateral signal has been found in multi-field models such as trapped inflation 
|44j . where there is no reduction in the speed of sound. In addition, it was shown recently 
that f^L may not depend explicitly on (? s in a weakly coupled extension of the EFT of 
inflation with reduced sound speed, depending instead on the Hubble scale and the UV 
cutoff of the theory [35]. Our result adds to the existing degeneracies as we have shown 
that even differentiating between canonical single field inflation models, with features in 
the potential, and non-canonical single field models may be more difficult than initially 
anticipated. However, differentiation is possible due to the tight constraints on the maximum 
amount of resonantly generated equilateral-type non-Gaussianity which arise from the non- 
observation of oscillating contributions to the 2-point function, and the structure of the 
resonant iV-point functions. In particular, this implies that detection of equilateral non- 
Gaussianity at a level greater than the PLANCK sensitivity of /jvl ~ 0(5) will rule out a 
resonant origin |20j. In case a detection of non-Gaussianity eludes PLANCK, future 21cm 
observations may provide a detection capability of Jnl ^ (9(0.01) down to the slow-roll 
level itself due to the extremely large number of modes (O(10 16 )) available [2T| 122] . For 
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equilateral-type non-Gaussianities of this magnitude there is then a degeneracy between a 
possible non-canonical and a resonant origin. 

Let us note in closing that we have consciously kept most of discussion at the level 
of effective field theory. Beyond the very general considerations laid out in Section 4 any 
concrete analysis of embedding quasi-equilateral non-Gaussianity from multiple resonant 
oscillatory contributions into string theory e.g. along the lines of axion monodromy has to 
await much fuller treatment in the future. 
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A Computing the mode functions and bispectrum 
A.l The mode functions 

We work in the comoving or unitary gauge where the inflationary perturbations are set to 
zero and spatial perturbations of the metric are written 



(A. 70) 



The Fourier transform lZ(k, t) can be written in terms of modes 



n(k, t) = n k (t)a(k) + nt(t)a\-k), 



(A.71) 
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-kr, for t the conformal time. For i < 1, TZ k (t) approaches a constant value lZ k ^ 



for which the mode function TZ k (t) satisfies the Mukhanov-Sasaki equation [121 H7], which 
for small e is given by [4"8] 

d 2 ft fc 2(1 + 2e + 5) dK k m 

-J-T-- L -r ± + n k = 0, (A.72) 

where x 

which is related to the primordial scalar power spectrum. In the slow-roll approximation 
the Mukhanov-Sasaki equation is easily solved to give the answer in (2.16) of [18J. But in 
our case, as in [IB], S/H is large compared to e and S and this SR approximation is invalid. 
Neglecting e and S and to leading order in e*, the Mukhanov-Sasaki equation becomes 

d 2 TZ k 2(1 + 5x{x)) dK k 



dx 2 



x 



dx 



o. 



(A. 73) 



for 



Si 



-3$>-sin 



+ Ci 



fi 



(A. 74) 



As in the single modulation case, the effect of Si is negligible at early times, for x ^> 



x tss a Vz, and 



where 1Z k °l is the value of 7Z k (x) outside the horizon in the absence of modulations and 
H^ji^x^j is a Hankel function. Similarly for x <C x res ^ i.e. at sufficiently late times, the effect 
of Si is again negligible and the solution must be of the form 



TZk(x) 



x^H^x) 



(A.75) 



TZh(x) 



K 



(o) 
k.O 



x ^ ( x ) 



2 k W 2 X3/ ' ^3/2 



As in ESI, 4+ } = 1 + 0((bl) 2 ) at late times, where c[ +) 



T c (+) 



(A. 76) 
1 + g(x) in the notation 



of [25]. (One can see this by solving (3.31) of [25] for each i) 
Thus we look for a solution of the form 



TZk(x) 



K 



(°) 
fc.O 



X ^ "^3/2 ( X ) 



c[ \x)i 



X ^^ "^3/2 ( X ) 



(A. 77) 



where cj^ (2;) vanishes at early times. Then the Mukhanov-Sasaki equation gives an equation 
governing the time evolution of c k (x) which is given by (for <C \f2e+ for all i) 

.Si(x) 



d 
dx 



-2i- 



x 



(A. 78) 
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This can be written as 



d 

dx 



-2ix 



d ( _) 
dx k 

d 

dx 



x) 
A{x) 



— ■ sm( 7; + 7; + — 



,/i J i 



x 



/i 



lnx) , 



lnx) 



(A.79) 
(A.80) 



Let A(x) = 'Y^ i A i {x), so that 



d 



! fc , Ci 



/i /i /i 



i J l r^k c i v - ' ; i 

cos( — + — H ; — in X) 



'2e 



ft fi ft 



is a solution. Then write 



2ix 



dx 



E4; 



(-), 



(A.81) 
(A.82) 

(A.83) 

<*. /,■/,■ /, -'■ (A ' 84) 

To find the mode functions outside the horizon, we use the expansions of the Hankel 
functions for small arguments: 



A(x) = J2Mx) = e~ 2tx ^ci-\x) = e- 



x) 



Then we find 



d_ 
dx 



c k,i ( x ) = 7==e cos(^- + — + -^^—Inx). 



4 /X\3/2 



2 x3/2 



2y / 7r \x J 
4 p\3/2 i /2\ 3/2 



(A.85) 
(A.86) 



We find 



K k {x) =R k % 



/i7T /" (f)k + Ci 

COS 



2^ 



ft 



*E 3fe * 



/iTT . f (pk + Ci 

: sin 



2^ 



(A.87) 



and 



(0)|2 _ |^(0) 1 2 



fc.Ol 



x/2^: 



1/2 



COS 



0fc + Cj 
/i 



(A.88) 
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A. 2 Mode functions in the bispectrum 



To see that we can use the unperturbed mode functions to calculate the bispectrum in Sec- 
tion 2.3 note first that at late times, for i < 1, the mode function in eq. (A. 88) approaches 
a constant, given by + cj^ ]. The contribution at early times, x 1, within the 

horizon, is 



n k (x) = n 



(0) 



-ixe lx +ic[ I 



x)xe~ 



-kfl 

where we have used the x>l limit of the Hankel functions 1181: 



-x 3/2 H^ 2 (x) « -zxe- ix 
and 7?-fc o is the unperturbed part of the mode function. We find 



^x^H^Jx) - 
2' 





~ ^-kfl Xe 


Ttk,l 




Ttkfl\ 







,(-) 

-fc 
d 

da; 



(-)/ 



(A.89) 



(A.90) 



(A.91) 
(A.92) 

(A.93) 



By eq. (2.18)) and eq. (A. 84) we see that this ratio is small for /j <C y2el (which is when 
the non-Gaussianity is appreciable). Thus the correction to the time derivative of the mode 
function is also small, and we can use the unperturbed mode functions and time derivatives 
in evaluating the three-point function. 



A. 3 Squeezed limit consistency relation 

In [TSJ the squeezed limit consistency relation found in [H] |23] was checked for resonant 
non-Gaussianity. It is straightforward to extend this to summed non-Gaussianities. The 
consistency condition can be phrased as [18] 

lim {K{ki,t)K{k*,t)K{k s ,t)) ~ -(27r) 3 5 3 (fc 1 + fc 2 + k*)\K™\ 2 \Tlf\ 2 d ^ } , (A.94) 
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where in the squeezed limit k\ ~ k 2 = k ^> k 3 . Furthermore, |7?.[ ' ) | 2 = 2?t 2 A^/A; 3 , witfj^ 



AUk) = AUK 
such that 



d \nA^(k) 
d Ink 



1/2 



sm 



^ k 
In 



neglecting slow roll corrections to n s = 1. Now we can write 



(A.95) 



(A.96) 



Um(7e(fc 1 ,t)7e(fc 2 ,t)7e(fc 3 ,t)) ~ (27r) 7 A^ — — 

A= 3 — >0 ^1^2% 



5 3 (fc! + fc 2 + fc 3 ) 



2fc 36V2tt /V2e* 
x — > — 
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sm 



In 
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(A.97) 



Comparing with eq. (2.25), to fulfill the consistency condition eq. (A. 94), Q res has to be given 

as 



Gres(k,k,k 3 ) 2k ^ 36*^2^ f y/2e 



k 2 k 3 



E 



fi 



1/2 



sm 



/ 2^ 1 *; 

In 



fi k* 



(A.98) 



in the squeezed limit. This is indeed the case up to a phase, as can be seen by taking the 



limit k 3 ->■ of eq. (2.32). 
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